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Topological quantum information processing relies on adiabatic braiding of nonabelian quasiparti¬ 
cles. Performing the braiding operations in finite time introduces transitions out of the ground-state 
manifold and deviations from the nonabelian Berry phase. We show that these errors can be elimi¬ 
nated by suitably designed counterdiabatic correction terms in the Hamiltonian. We implement the 
resulting shortcuts to adiabaticity for simple protocols of nonabelian braiding and show that the 
error suppression can be substantial even for approximate realizations of the counterdiabatic terms. 


Introduction .—It is envisaged that the information 
processing of topological quantum computers relies on 
adiabatic braiding of nonabelian quasiparticles [IHS] ■ Ex¬ 
changing two nonabelions does not leave the quantum 
state unchanged, possibly up to a sign (as for fermions 
or bosons) or phase (as for abelian anyons) factor, but 
rather effects a unitary rotation in a degenerate subspace 
of ground states. The ground-state degeneracy grows 
exponentially with the number of nonabelian quasipar¬ 
ticles, and quantum information processing corresponds 
to manipulating the system’s ground state by braiding 
operations. Majorana bound states in topological su¬ 
perconducting phases constitute the simplest example of 
such nonabelions [B] , and there has been considerable ex¬ 
perimental effort towards realizing a possible hardware 
[THIS], following a series of theoretical proposals [T4H22] . 

Topological quantum information processing is im¬ 
mune to local sources of decoherence when braiding is 
performed adiabatically [T]. Quite generally, adiabatic¬ 
ity is protected by the gap of the underlying topological 
phase. Here we want to ask the question whether it is 
possible to realize the exact adiabatic quantum dynamics 
of the braiding operation, albeit in a finite time interval. 
There are obvious motivations why this would be desir¬ 
able: First, any topological quantum computer would 
operate at a finite clock speed which necessarily entails 
possibly small, but nonzero errors. Second, a topolog¬ 
ical quantum computer would presumably have to op¬ 
erate faster than parasitic decoherence processes such as 
quasiparticle poisoning or deviations from perfect ground 
state degeneracy originating in the finite spatial extent 
of the Majorana quasiparticles. In both cases, such a 
scheme could then be used to offset the incurred errors - 
enabling longer computations or higher clock speeds. 

Demirplak and Rice [23] as well as Berry [24] intro¬ 
duced a protocol that emulates the adiabatic dynamics of 
any nondegenerate Hamiltonian Hq (t) as the exact quan¬ 
tum dynamics in finite time. This scheme is known al¬ 
ternately as transitionless quantum driving or shortcut 
to adiabaticity. The prize that comes with the short¬ 
cut is that the adiabatic quantum dynamics of Hgit) is 
generated by a Hamiltonian H(t ), which differs from Hq 


by counterdiabatic terms Hi{t). This shortcut to adia¬ 
baticity does not apply directly to the adiabatic braid¬ 
ing of nonabelian quasiparticles because of the associated 
ground-state degeneracy. Here, we first generalize this 
scheme to systems with degenerate manifolds of states 
where adiabatic dynamics generates nonabelian Berry 
phases. Then, we apply this generalized shortcut to non¬ 
abelian statistics, using a simple model for braiding of 
Majorana bound states. Within this model, the braiding 
of Majorana zero modes is based on judiciously chosen 
temporal variations of the couplings between a number 
of Majorana end states. We find that shortcuts to non¬ 
abelian braiding can be implemented by introducing a 
small number of additional local couplings. 

Shortcuts to adiabaticity for degenerate systems. —The 
exact quantum dynamics of a Hamiltonian H(t) is gener¬ 
ated by the corresponding time-evolution operator U(t) 
which satisfies the Schrodinger equation 

idtlfit) = H{t)U{t). (1) 

Thus, we can give an explicit expression for the Hamilto¬ 

nian H (t) generating any prescribed quantum dynamics 
U{t), 

H{t) = i[dtU]lil (2) 

The shortcut to adiabaticity [23122] follows by inserting 
into this expression the adiabatic time-evolution operator 

U{t) = ^e-*/o (3) 

n 

for the Hamiltonian Ho{t), with instantaneous eigenval¬ 
ues En{t), instantaneous eigenstates |■^/'n(I)), and Berry 
phase 7 n(t) = i dt' One finds that 
H{t) = HQ{t) + Hi{t) with the so-called counterdiabatic 
terms [23!^ 

Hl{t) = i'^(\dtfin){lfn\ - |V’n)('0n|l9tV’n)(^/’n|) • (4) 

n 

Such shortcuts to adiabaticity have recently been im¬ 
plemented experimentally for effective two-level systems 
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arising in trapped Bose-Einstein condensates [25] and for 
the electron spin of a single nitrogen vacancy center |26j . 

Following Wilczek and Zee we now consider a 

Hamiltonian 11^(1) whose instantaneous spectrum de¬ 
fined through 

Mmiit)) = (5) 

includes one or more sets of states |'i/' 2 (^)) which remain 
degenerate for all t. Here, a = 1,... ,dn labels the states 
within the degenerate subspace n of multiplicity dn- 
We first define the adiabatic solution of the 

time-dependent Schrodinger equation 

idMlit)) = HQ{t)\ri^{t)) (6) 

with initial condition |?? 2 ( 0 )) = l^a(O))- the adiabatic 
limit, the time-evolved state need not remain parallel to 
IV'a(t)) but will in general be a linear combination of all 
basis states within the degenerate subspace, 

\rim) = Y.^2p\m))- ( 7 ) 

d 


Inserting this expansion into the time-dependent 
Schrodinger equation yields an equation for the coeffi¬ 
cient matrices C/", 


idtU^ = U^{A^-Er,l), (8) 


where denotes the nonabelian Berry 

connection m- This is solved by 

^ g-ifo 

in terms of time-ordered exponentials. 

The adiabatic time evolution of the Hamiltonian Hq (t) 
follows the time-evolution operator 


m = E = E 

n,a n,a^ 

( 10 ) 

Now we use Eq. ([^ to derive the Hamiltonian H{t) for 
which this is the exact time-evolution operator. Inserting 
Eq. (10) into we obtain the shortcut to adiabaticity 
(all quantities evaluated at time t) 


H = 


n 

+[(c/”)^c/"]/3a|5tC)(V'^l} (11) 


The second term in H simplifies due to unitarity of 
I/", ([/”)ll/” = 1. Combining unitarity and Eq. (j^, 
we also have i(t7")lt/" = (E^l — H") which simplifies 
the first term. With these identities, we readily find 
H{t) = Ho{t) + Hi{t) with 


iJi = iE 

n 


Eiw(d-Ei^2)(ci5t^?)(V'^i ■ 

a 0.13 

( 12 ) 


These counterdiabatic terms generalize the shortcut to 
adiabaticity to systems with degenerate spectra and non¬ 
abelian Berry connections. 

Majorana systems. —In view of topological quantum 
information processing, we specifically consider the coun¬ 
terdiabatic terms for a Bogoliubov-de Gennes Hamilto¬ 
nian in Majorana representation, 

Eto ~ ^ ^ ^ ^n'yn,2a — l^n,2a- ( 1 ^) 

no 


Here, both e„ and the jn.a are explicitly time dependent 
and associated with the instantaneous Hamiltonian. The 
instantaneous many-body spectrum of Hq contains de¬ 
generacies whenever an eigenenergy e„ vanishes or when 
one or several nonzero e„ are degenerate. The Majorana 
eigenmodes associated with e„ are denoted by 7 „_q, where 
a takes on 2N values for an TV-fold degenerate energy e^. 
The counterdiabatic terms Hi guarantee that the time 
evolution generated by the full Hamiltonian Hq+Hi does 
not take the Majorana eigenmodes jn.a out of the sub¬ 
space n. At the same time. Hi should not alter the time 
evolution within these subspaces. In the supplementary 
material m we show that these conditions yield 

Hi = ~ ) '7n,a7n,a ~ g ^ ( 7ra,a j 7n,/3 }7n,/3 ■ (14) 

no n,o 0 


This result complements the counterdiabatic terms in 
first quantization in Eq. (12). 

Application to nonabelian braiding. —A minimal model 
for nonabelian braiding starts from a Y-junction of 
three one-dimensional topological superconductors, la¬ 
beled wire I, 2, and 3 [SHUfiCi] . as illustrated in Fig. [^a). 
If all three arms are in the topological phase, there are 
four Majorana bound states in this system. Three of 
these are located at the outer ends of the three wires, 
with Bogoliubov operators labeled "fj for wire j, and a 
fourth Majorana mode 70 is located at the junction of 
the three wires. As long as the three arms have a finite 
length, these outer Majorana bound states hybridize with 
the central Majorana and the system is described by the 
Hamiltonian 
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= * E ^“To7a (15) 

a=l 


This Hamiltonian couples the central Majorana 70 to a 
linear combination of the outer three Majoranas. We can 
thus readily bring it to the form of Eq. (13), 

Hq = z/ia7o7a, (16) 


with 7A = (1//ja) Z)a=i ^a7a and Ha = [Af + + 

For any choice of the couplings A^-, there are also 
two linearly independent combinations of the outer Ma¬ 
joranas which do not appear in the Hamiltonian and thus 
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FIG. 1. (a) Y-junction with central Majorana 70 and three 
outer Majoranas 7 j {j = 1,2,3). The outer Majoranas are 
coupled to the inner Majoranas with strength Aj. (b) Basic 
step of the braiding procedure, moving a zero-energy Majo¬ 
rana from the end of wire 1 to the end of wire 3 by varying 
the Aj. Dark (light) wires indicate zero (nonzero) couplings 
Aj. Dark red circles correspond to zero-energy Majoranas, 
light blue circles indicate Majoranas acquiring a finite energy 
by coupling. In the intermediate step, the zero-energy Majo¬ 
rana is delocalized over the three Majoranas along the light 
wires, (c) Three steps as described in (b) result in braiding 
the zero-energy Majoranas 71 and 72 . 


remain true zero-energy Majoranas. Due to these zero- 
energy modes, the two eigenvalues of Hq are each doubly 
degenerate. Specifically, when just one of the couplings 
Aj is nonzero, these two zero-energy Majoranas can be 
identified with the Majoranas located at the ends of those 
wires with zero coupling. 

We assume that we can change the couplings Aj as 
a function of time. We can now imagine the following 
braiding procedure [IHUSO]- Initially, only A 3 is nonzero. 
Then, 71 and 72 are zero-energy Majoranas. In a first 
step, we move a zero-energy Majorana from the end of 
wire 1 to the end of wire 3, without involving the zero- 
energy Majorana 72 as shown in Fig. [^b). To this end, 
first increase Ai to a finite value. The zero-energy Ma¬ 
jorana originally located at the end of wire 1 is now delo¬ 
calized and a linear combination of 71 and 73 . We then 
localize the Majorana zero mode at the end of wire 3 by 
reducing A 3 down to zero, leaving only Ai nonzero. The 
braiding process is completed by two completely analo¬ 
gous moves (see Fig.j^c)): We first move the zero-energy 
Majorana from the end of wire 2 to the end of wire 1, and 
finally the zero-energy Majorana from wire 3 to wire 2. 
The combined effect of this procedure is to exchange the 
initial zero-energy Majoranas at the ends of wires 1 and 
2. One can check easily [30] that the change of the state 
of the system under this adiabatic exchange is described 
by the familiar braiding matrix C /12 = exp(i 7 r 7 i 72 / 4 ). 


When performing this exchange operation over a finite 
time interval, there will be corrections to the adiabatic 
time evolution. We can now apply one of the nonabelian 
shortcut formulas in Eqs. (12) or (141. As shown in the 


(a) 





FIG. 2. (a) Minimal implementation required for braiding 
with shortcut protocol. The additional couplings needed 
for the shortcut protocol are shown in blue, (b) Wire net¬ 
work with many Majoranas allowing for pairwise exchanges 
of neighboring Majoranas including shortcut protocol. Im¬ 
plementing the shortcut merely requires the addition of local 
couplings within the network. 


supplementary material EH, we obtain 

Hi = ^7a7a = ^ ^(A^A„ - A„A^) 7 „ 7 ^. (17) 

^ a</3 

Thus, the shortcut is based on additional couplings be¬ 
tween the outer Majoranas, while the adiabatic braiding 
protocol only uses couplings between the central and the 
outer Majoranas, see Fig. [^a). Specifically, during the 
basic step of moving a zero-energy Majorana from the 
end of wire i to wire j, only the couplings A^ and Aj 
are nonzero. According to Eq. ( |17| ), performing this step 
accurately in finite time merely requires the additional 
coupling between 7 ^ and 7 ^. 

Practical implementation .—There has been consider¬ 
able work on how to implement braiding based on one¬ 
dimensional superconducting phases [28H301132H34] . The 
couplings of the Majoranas can, e.g., be varied by chang¬ 
ing the length of the intervening topological section. 
However, this may not be easily compatible with the ge¬ 
ometric constraints imposed by the shortcut protocol, cf. 
Eig. [^a). A better approach may be to vary the mag¬ 
nitude of the topological gap. Both methods control the 
overlap of the Majorana end states and hence their cou¬ 
pling. Physically, this can be achieved, say in quantum- 
wire based realizations, by changing the chemical poten¬ 
tial by means of a gate electrode [28] or a supercurrent 
in the adjacent s-wave superconductor [35] . 

More controlled variations of the Majorana couplings 
may be possible by exploiting charging effects [30] or by 
quantum dots |32| . Eor simplicity, assume that the quan¬ 
tum dot has a single level which is tunnel coupled to the 
ends of two topological wires with their Majorana end 
states. When the dot level is far from the Fermi energy, 
there is essentially no coupling between the adjacent Ma¬ 
joranas. Conversely, when the dot level is close to the 
Fermi energy, the Majoranas become strongly coupled. 
This approach modifies the coupling of the Majoranas 
by conventional gate control of a quantum-dot level and 
is also compatible with the geometric constraints of the 
shortcut protocol. 
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So far, we have focused on the exchange of two Ma- 
joranas within the minimal setting of a Y-junction. Of 
course, one can readily imagine a scheme in which there 
is an entire keyboard of Majoranas and any two neigh¬ 
boring Majoranas can be readily braided. Importantly, 
amending this scheme to implement the counterdiabatic 
terms merely requires additional local couplings as shown 
in Fig.j^b). 

Robustness. —The manipulation of the quantum state 
is independent of the precise braiding path as long as 
the exchange is performed adiabatically. In contrast, 
the diabatic corrections are sensitive to the details of 
the braiding protocol. Consequently, the counterdiabatic 
terms 0 are not topologically protected, depend on the 
specifics of the braiding path, and for full effect, have to 
be implemented exactly for a given Ho{t). 

However, we find that one can reach substantial re¬ 
ductions in the diabatic errors even when the shortcut 
protocol is implemented only with reasonable accuracy. 
We have computed the diabatic errors numerically, both 
for the bare braiding protocols and for approximate im¬ 
plementations of the counterdiabatic terms. Specifically, 
we consider the diabatic errors for 

Hxit)=Ho{t) + XHi{t). (18) 

For A = 1, the counterdiabatic terms exactly compensate 
the diabatic corrections for any duration of the braid¬ 
ing protocol. As approximate implementations of the 
counterdiabatic terms, we consider relative errors of 10 % 
(A = 0.9) and 30% (A = 0.7). We compute both the tran¬ 
sition probability out of the degenerate subspace and the 
accumulated deviation from the adiabatic Berry phase. 

Implementing the basic step [shown in Fig. [^b)] of 
the braiding protocol in Fig. [^c) by Ai(t) = Asin(/?(f) 
and A 3 (t) = Acos with (f{t) increasing from 0 to 
7r/2, both the transition probability and the phase error 
exhibit a power-law dependence on the protocol duration 
T. The power law depends on the specific choice for ^p{t). 
Choosing the latter such that the derivative vanishes at 
the end points yields a T~^ dependence. Corresponding 
numerical results are included with the supplementary 
material [3T] . Interestingly, we find similar results for the 
protocol given in Ref. m, in which one initially increases 
Ai, leaving A 3 constant, and then reduces A 3 to zero in 
a second step [3T] . 

Exponentially small transition rates can be realized by 
choosing Ai(t) = Asin^(/?(t) and A^^t) = Acos^ 

Now the gap assumes a minimum during the protocol as 
in the familiar Landau-Zener process. For the numerical 
calculation presented in Fig. we have chosen ip{t) to 
have a smooth derivative. The diabatic transition rate 
is indeed exponential in the protocol duration which is 
somewhat conterintuitive as the transition rate actually 
decreases relative to the previously discussed protocols 
although the gap is smaller. The phase error also exhibits 
exponential scaling as shown in the inset of Fig. 



duration (1/A) 


FIG. 3. Diabatic errors vs duration of braiding protocol for 
the transition probability out of the degenerate subspace of 
the initial state. The inset shows the phase error relative 
to the nonabelian Berry phase. For both quantities, curves 
are shown in the absence of counterdiabatic terms [A = 0 in 
Eq. (18 1 ] and with counterdiabatic terms with 10% (A = 0.9) 
and 30% (A = 0.7) relative error. There would be no diabatic 
error if the counterdiabatic errors were implemented exactly. 


An exact implementation of the counterdiabatic terms 
fully corrects for these errors. As can be seen from Fig. 
a suppression by two orders of magnitude merely requires 
an implementation which is accurate at the 10 % level. 
Even a very rough implementation at the 30% level still 
substantially reduces the errors. More generally, we find 
that the relative error scales approximately as (1 — A)^ 
with the accuracy of the implementation of Hi. It is 
also worth noting that the approximate counterdiabatic 
terms suppress the diabatic error, but do not modify its 
scaling with protocol duration. 

Conclusions. —In summary, we have generalized the 
concept of shortcuts to adiabaticity to nonabelian Berry 
phases and showed how this can in principle be used to 
implement nonabelian braiding operations exactly in a 
finite time. Such protocols can substantially improve the 
accuracy of braiding operations performed in a finite time 
interval. It is interesting to note that our scheme bears 
some resemblance with the concept of quasi-adiabatic 
continuity for topological phases [36j. 

In this work we have focused on a simple model of non¬ 
abelian braiding which excludes the quasiparticle contin¬ 
uum. The current protocols are therefore useful when¬ 
ever there is a separation of scales between the finite- 
energy subgap states and the magnitude of the topolog¬ 
ical gap. Including the quasiparticle continuum is an 
interesting problem for future research. It should also 
be interesting to extend the current considerations for 
Majorana zero modes to more exotic nonabelian quasi¬ 
particles. 
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SUPPLEMENTARY MATERIAL 


Derivation of Eq. (14) in the main text 


Here we derive general expressions for the counterdiabatic correction for non-interacting systems in terms of the 
(instantaneous) eigenmodes ')n,a of the original (time-dependent) Hamiltonian in a Majorana representation, 


Hq — i £n7n,2a-l'ln,2a- 


(19) 


Note that both the eigen-Majoranas 7 „_q and the eigenvalues e„ are time dependent. Degeneracies of the many-body 
spectrum can arise when one or more of the eigenvalues e„ vanish or when some nonzero e„ is degenerate, independent 
of time. The various Majorana operators associated with each single-particle eigenvalue e„ are labeled by the index 
a. If the single-particle eigenvalue e„ is iV-fold degenerate, a takes on 2N different values, a = 1,..., 2N. 

A direct derivation of the counterdiabatic terms based on the general Eq. (12) in the main text is cumbersome. Here, 
we choose to proceed as follows. The counterdiabatic terms suppress transitions out of the degenerate subspace but 
leave the dynamics within the degenerate subspace as governed by the nonabelian Berry connection of the original time 
evolution unchanged. Thus, we can determine the counterdiabatic terms Hi uniquely from the following constraints: 

(a) Hi has no matrix elements which act within the degenerate eigenspaces of Hq. This ensures that Hi affects 
only transitions between states with different energies. 

(b) The time evolution of the 'yn,a with respect to the full shortcut Hamiltonian H — Hq + Hi does not include any 
transitions between different degenerate subspaces. 

To implement these constraints, we note that the time evolution of the Majorana operators is governed by the 
Heisenberg equation of motion 


d'yn,a _ d"fn,o 
dt 


dt 


+ i[Ho + Hi, 7n,Q.], 


( 20 ) 


where the first term on the right-hand side accounts for the explicit time dependence of the Majorana operators. This 
term allows for the expansion 


9jn,a _ 

~ / . ^nmlrn,p- 


( 21 ) 


m,/9 


Here, the coefficients can in principle be expressed in terms of the instantaneous eigenfunctions and their time 
derivatives. It turns out, however, that we do not need these explicit expressions for the present purpose. To proceed, 
we also write the counterdiabatic terms in the general form 


Hi — i ^ ^ ^ ' ^n,’m7n,a7m,/3- 


( 22 ) 


n,m a,P 


Thus, it is our goal to derive the coefficients which satisfy the antisymmetry relation 
Implementing the contraints (a) and (b), we demand that the 'yn,a satisfy the time evolution 


d'Jn 


dt 


= +i[Ho,Jn,a\- 


(23) 


Here, the first term of the right-hand side contains only those terms of 


terms in 


d-Jn 


^7n,c 

dt 


that belong to the same subspace n. All 


dt 


which belong to different subspaces must be cancelled by the counterdiabatic terms Hi. 


The desired shortcut time evolution in Eq. (23) satisfies 

d^n, 


dt 


Tm,/3 


= 0 


(24) 


for m ^ n. Inserting the Heisenberg equation of motion (201 into this condition, we obtain 


Ah^P _ Aj^pa I 


d^n,c 


) ^m,P r — 0 


(25) 
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for m ^ n. Using the antisymmetry property of the yields 


'^> = 4 EE 


n,m 

(n^m) 




5 Tm,/3 f • 


Finally, we write this as 


Tj _ J 1 i / djn 

-Hi — 7m,^ 7m,/3, ^ 7n,a 7ri,/3 7n,/3; f 7n,c 


n,m oi^ 


and use the relation 


Y^Y^ I d^n,a\ _ ^d^n, 

/ , / , 7m,/3 ^ 7m, - ? — 2 


m /3 


dt 


dt 


(26) 


(27) 


(28) 


to obtain Eq. (14) of the main text. 

In the following, we derive the counterdiabatic terms for the braiding procedure given in Eq. (17) of the main text 
using this general result as well as a more basic approach starting with Eq. (12). 


Derivation of Eq. (17) using the general Majorana counterdiabatic terms 

Here, we derive Eq. (17) using Hi in the Majorana operator representation as given in Eq. (14). The braiding 
Hamiltonian in Eq. (16) is 

Ho =ih A (i)7o7A (t) . (29) 

The system comprises a mode of energy hA{t) associated with the two Majorana operators 70 and 7 a( 0 as defined in 


the main text. In addition there are two Majoranas 7 Ait) and 73 ( 1 ) which remain uncoupled by Eq. (29) and form a 
zero-energy mode. Using the identity 


{7a, 7 / 3 } = -{7/3, 7a} 


(30) 


(with the shorthand 7 = ^) we can write the time derivatives (suppressing time arguments) in the most generic form 
as 


7A = Vi7b +»?27a, 

7 b = -r]i7A + r737A, 

7A = -m7A - m7B- 

with real coefficients rji. The first term in Hi can be written as 

F 7n,a7n,a = 7a7a + 7b7b + 7a7a = 27A7A + ‘2r]i7B7A 


(31) 


(32) 


The second contribution to Hi subtracts all terms within a degenerate subspace and thus eliminates the term ~ 7b7a- 
Thus we obtain 


Hi = ^7A7A = ^ ^a A^7„7/3 


(33) 


A 


a/3 


as given in Eq. (17). 


Derivation of Eq. (17) using the spin construction 


We can alternatively derive Eq. (17) using the general formulation of the counterdiabatic terms in Eq. (12). To this 
end, we introduce conventional fermionic operators through 

Cl = ^(71-*72) ; C2 = ^(70-173)- 


(34) 








Using the inverse relations 


7i = Cl + c| ; 72 = i{ci - c|) ; 73 = i{c 2 - cl) ; 70 = C 2 + cl, 


(35) 


we can write Hq in terms of ci and C 2 
3 

Ho = iY^ Aj 7o7j 


= iAi(c 2 Ci - c|c2 + C2C1 - c\cl) - A 2 {clci + c|c2 + C2C1 + c\cl) - A 3 ( 2 c 2 C 2 - 1). 
Specifically, we write the Hamiltonian in the basis {|00), |11), |10), |01)}, where the basis states are defined as 


(36) 


111) = cjcJlOO) , |10) = cjlOO) , |01) = cj|00) 


with ci|00) = C 2 IOO) = 0. This yields 


/ A 3 
—iAi — A 2 
0 

V 0 


iAi — A 2 
—A 3 
0 
0 


0 0 \ 

0 0 

A3 —lAi — A2 

iAi — A2 —A3 J 


(37) 


(38) 


The block-diagonal structure originates from the conservation of fermion-number parity. In fact, it is easy to show 
that the Hamiltonian H commutes with the parity operator 


P = 7o7i7273- 


(39) 


The top-left block i7even = A^Tz — AiTy — A 2 TX corresponds to even fermion parity, while the bottom-right block 
Hodd = A^Tz + AiTy — A 2 TX has odd fermion parity. Here we have defined Pauli matrices Ti within the even and odd 
subspaces. If we also define Pauli matrices tTj in the even-odd subspace, then we can write 

Ho = AoTz - AiTyTTz - A 2 TX (40) 

for the overall Hamiltonian H. Expressing iJgven and H^dd in terms of Pauli matrices makes it obvious that 
these Hamiltonians take the form of a spin Hamiltonian in magnetic fields Seven = (—A 2 ,—Ai,A 3 ) and Sodd = 
(—A 2 , Ai, A 3 ), respectively. The degeneracy due to the presence of the Majorana modes implies that the two sub¬ 
spaces have the same eigenvalues. At the same time, the spectrum for each subspace by itself is non-degenerate. 

In order to evaluate the counterdiabatic terms, it is useful to eliminate the time derivatives of the states from Eq. 
(12). To achieve this, we first multiply the first term on the right-hand side of Eq. (12) by 1 = S /3 

from the left and obtain 


m^n oc^ 


(41) 


Taking the time derivative of Eq. (5) and multiplying from the left by ('0™|, one finds 

{^f\dtHo\ro.) 




E — E 

-^n 


for n m. Inserting this into Eq. (41) yields 


w;i- 


m^n 0.(3 




(42) 


(43) 


Using this expression, the counterdiabatic terms Hi can be conveniently derived. 

To do so, we temporarily perform rotations within the even and odd subspaces such that Hq {t) in Eq. (401 involves 
only the Tz term. Then, the eigenstates in the even and odd subspaces are simply the “spin-up” and the “spin-down” 
states. Using 


dtHo = AsTz - AlTyTTz - A2TX, 


(44) 
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duration (1/A) 

FIG. 4. Diabatic errors vs duration of the braiding protocol defined by Ai(f) = Asin^(i) and A 3 (f) = Acos(^(f) for the 
transition probability out of the degenerate subspace of the initial state. The inset shows the phase error relative to the 
nonabelian Berry phase. The excitation gap remains unchanged during the entire braiding protocol, which results in a power- 
law dependence on the duration. The function is chosen to have zero derivative at both endpoints. For both transition 
probability and phase error, curves are shown in the absence of counterdiabatic terms and with counterdiabatic terms with 
10% and 30% relative error. There would be no diabatic error if the counterdiabatic errors were implemented exactly. 


we then find that 


2(Af + Ai + Ai) 


|(A3Ai — AiA3)Ta;7r2 -h (A2A3 — A3A2)tj, — (A 1A2 — A2Ai)rz7r2| . 


This can be readily expressed in terms of the original Majorana operators. Indeed, we have the identities 

*Ai27i72 = -Ai 2 ( 2 c|ci - 1 ) 

*Ai37i73 = Ai3(c2Ci -I- c\cl - cjci - c\c 2 ) 

*A237273 = *A23(c2Ci - c\cl + c\c2 - 4ci), 


or, in the basis specified above. 


iAi27i72 = Ai2Tz7rz 
iAi37i73 = Aist^jTTz 

iA237273 = -A23Ty. 


Thus, we finally find 


Hi = 


2 (A? + Ai + Ai) 

in terms of the original Majorana operators. 


|(A2Ai — AiA2)7i72 + (A3A1 — AiA3)7i73 -I- (A3A2 — A2A3)7273| 


(45) 


(46) 


(47) 


(48) 


Numerical calculation of the robustness 

In this section we provide details of the numerical calculations. For completeness we also include numerical results 
for the transition probability and Berry phase errors of the non-exponential protocols mentioned in the main text. 
Due to the conservation of fermion-number parity the nonabelian Berry phase takes the form exp(z 7 rz), where is a 
Pauli matrix in parity space. Performing the braiding protocol adiabatically yields 7 = 7 r/ 4 . For finite durations we 
numerically compute the Berry phase as 7 = arg[(4'e(T)4'e(0))/(4'o(T)4'o(0))]/2, where 4'e/o(0 denotes the ground 
state wavefunction at time t with even (odd) parity. The diabatic phase error is I 7 — 7 r/ 4 |. 

We first consider the protocol with the basic step Ai(t) = Asinip{t) and A 3 (t) = Acos(p{t). When has zero 
derivative at both endpoints ((/? = 0 for t = 0 and ip = 7 r /2 for t = T/3), the transition probability scales as 
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FIG. 5. (a) Couplings Aj and gap /ia during the braiding process for the protocol in Ref. [30]. (b) Diabatic errors vs duration 
of braiding protocol. Similar to the protocol in Fig. the error has a power-law dependence on the duration. The derivative 
of jumps at the end points and therefore the error scales as T~^. 


with the protocol duration T. This is shown in Fig. Specifically, we have chosen (p(t) = ( 7 r/ 2 )[ 3 ( 3 t/T)^ — 2{3t/T)^] 
between 0 and t = T/3 for this calculation (as well as for the one in the main text). When the derivative jumps at 
either (or both) end points (as for the simplest choice (p{t) = 3TTt/2T), we find the errors to decay even more slowly, 
namely as T“^. Note that both the transition probability out of the degenerate subspace and the phase error of the 
topological qubit scale in the same manner with T, see inset of Fig. 

Interestingly, the same dependences are found for the protocol given in Ref. [30] and displayed in Fig. [^a). The 
initial step of the braiding operation is effected by increasing Ai first at constant A 3 . The latter is reduced to zero 
only subsequently. In this protocol, the gap increases and takes on a maximum halfway through this basic step. 
Nevertheless, the diabatic errors still vary as a power law of T. Fig.j^b) shows corresponding numerical results. Here 
we chose a linear protocol, ip(t) = 37 rt/ 2 T, in which the derivatives of ip{t) do not vanish at the end points. 































































